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ABSTRACT 



Protons of energies up to ~ 10 2U eV can be subject to significant deflection and 
energy dependent time delay in lage scale extragalactic or halo magnetic fields of 
strengths comparable to current upper limits. By performing 3-dimensional Monte 
Carlo simulations of nucleon propagation, we show how observations of arrival direction 
and time distributions can be used to measure the structure and strength of large-scale 
magnetic fields, and constrain the nature of the source of ultra-high energy cosmic 
rays. 



Subject headings: cosmic rays - magnetic fields 



- 2 - 



1. Introduction 

If extragalactic magnetic fields (EGMFs) exist, they are at a level below detectability with 
presently available techniques. Faraday rotation measures of distant powerful radio-sources give 
upper limits to extragalactic fields of B rms y/T c £ 10~ 9 G Mpc 1//2 , where l c denotes the reversal 
length of the field (see, e.g., Kronberg 1994). However, EGMFs below these current limits are 
of significant interest to cosmology, galaxy and star formation, and galactic dynamos (see, e.g., 
Kronberg 1994; Olinto 1997). 

The propagation of UHECRs is affected by the presence of EGMFs of strength 
1(T 12 G £ B ims £ 1(T 9 G and/or by Galactic halo fields 1(T 8 G £ B H ^ 10~ 6 G. Protogalactic 
fields in the former range are actually expected if the galactic magnetic field cannot be explained 
by a galactic dynamo (Kulsrud & Anderson 1992). In this letter, we study how EGMFs and 
the halo magnetic field, which we collectively name large-scale magnetic fields (LSMFs), affect 
UHECRs with energies E ^ lOEeV (EeV = 10 18 eV). We simulate the propagation of ultra-high 
energy nucleons in the LSMF and the cosmic microwave background (CMB), and show how the 
resulting angle-time-energy images of these UHECRs, could be used, in turn, to probe LSMFs 
with strengths in the above range. Suitable UHECR statistics could be achieved with future 
experiments, such as the Japanese Telescope Array (Teshima et al. 1992), the High Resolution 
Fly's Eye (Al-Seady et al. 1996), and the Pierre Auger Project (Cronin 1992), that have the 
potential to detect ~ 10 — 10 3 particles per source with E ^ 10 19 eV, over a ~ 5yr period. 
Although the origin and nature of UHECRs are not known, they are likely to be generated in 
extragalactic sources {e.g., Bird et al. 1995; Yoshida et al. 1995; Hayashida et al. 1996), such as 
powerful radio-galaxies {e.g., Rachen & Biermann 1993), cosmological 7— ray bursts (Vietri 1995; 
Waxman 1995; Milgrom h Usov 1996), and/or topological defects {e.g., Bhattacharjee, Hill & 
Schramm 1992; Sigl 1996). Our results can also be used to discriminate between models of the 
origin of UHECRs. 

Previous studies of the effect of the LSMF on UHECRs have included a discussion of the 
energy dependent deflection and time delay for extragalactic UHECRs (Cronin 1992; Waxman & 
Miralda-Escude 1996; Medina Tanco et al. 1997). In addition, the effect of EGMFs on secondary 
7-rays produced in the interaction of ultra-high energy protons with CMB photons can also probe 
EGMFs below current upper limits: Plaga (1995) and Waxman & Coppi (1996) proposed to 
use the arrival time delay of secondary 7-rays in the TeV range to probe EGMFs of strengths 
B ims £ KT 15 G, while Lee, Olinto & Sigl (1996) showed how fields of strength B ims ~ 1(T 9 G affect 
the 7— ray spectrum around 10 EeV due to electron synchrotron losses. 

2. Angle-time-energy images 

Nucleons propagating in intergalactic space with energies ^ 10 EeV are mainly subject to 
scattering on the LSMFs and pair production on the CMB (for protons), as well as photopion 
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production on the CMB. Pair production dominates the energy loss below ~ 70EeV which we 
include as a continuous loss process (Chodorowski, Zdziarski, & Sikora 1992). Above ~ 70EeV 
photopion production dominates and gives rise to the Greisen-Zatsepin-Kuz'min cutoff (hereafter 
GZK cutoff; Greisen 1966; Zatsepin & Kuzmin 1966). We model photopion production as a 
stochastic process where the chance of interaction is drawn at random, as described in Lee (1996), 
and Sigl, Lee, & Coppi (1996). We model the EGMF as a gaussian random field, with zero mean 
in Fourier space, and power spectrum: (B 2 (k)) oc k riB for k < 2ir/l c , and (B 2 (k)) = otherwise. 
Here, l c characterizes the cutoff scale of magnetic fluctuations; we choose l c = 1 Mpc as a fiducial 
value for the EGMF. The power spectrum is normalized via: B 2 ms = V/(2t:) 3 J d 3 kB 2 (k). We 
leave the detail description of our Monte-Carlo code to Sigl, Lemoine, & Olinto (1997). 

A simpler version of the present study was carried out by Medina Tanco et al. (1997). These 
authors described the field as an assembly of randomly oriented bubbles of constant field, with 
diameter equal to the coherence length, and neglected the stochastic nature of pion production. 
Finally, they considered the limit DOe 3> l c , where D denotes the distance to the source, and 
9e the r.m.s. deflection angle at energy E; this limit is obtained only for extreme choices of the 
magnetic field strength and configuration for E ;> 10 20 eV (see below). The limit DOe <C l c is not 
only the most probable for E > 10 19 eV, it also is the most difficult to treat (see Sigl, Lemoine, & 
Olinto (1997). The importance of distinguishing the limits DOe 3> l c and DOe <C l c , was pointed 
out by Waxman & Miralda-Escude (1996). In the latter limit, all nucleons emitted by the source 
and captured on the detector have essentially experienced the same magnetic structure along 
their paths, hence the scatter around the mean of the correlations between the time delay te, the 
deflection angle Oe, and the energy E, is of order 1%. In the opposite limit, D9e 3> l c , the relative 
widths of these correlations are of order unity. 

Our simulations generate data consisting of the energy, arrival time, and angular direction 
for each UHECR emitted from the source. For general reference, we can estimate the average 
deflection angle Oe and the average time delay te — D0 E /2c induced on a proton of energy 
E <; 50EeV over a distance D in the limit DOe > l c - Using (Waxman & Miralda-Escude 1996): 

TE " 2 - 5 (ItS) (iom^c) (ioeW) {w^g) (iM^) yr ' (1) 

and E ^ 0.014° (D/10 Mpc)- 1 / 2 ^/! yr) 1/2 . The Faraday rotation bound B rms ^/T c < 
lO^GMpc 1 / 2 , implies t e ^ 3 x 10 4 (L>/10Mpc) 2 (^/10EeV)" 2 yr, and E <, 2.5°(£>/10 Mpc) 1 / 2 
(.E/10 EeV) -1 , assuming ub — 0. In what follows, we address the possible observables for different 
cases. In §2.1 and §2.2, we discuss the case of an extra-galactic bursting source. In §2.3, we discuss 
the case of a continuously emitting source. In these sections, we consider the case in which the 
time delay induced by the extra-galactic magnetic field dominates over that of the halo magnetic 
field of our Galaxy. In §2.4, we discuss how the domination by a halo field would modify our 
conclusions. We use the notations tiqo = TE=woEeV, #ioo = #E=ioo.EeV- 
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2.1. Observable time delays 

The best case for probing large scale magnetic fields as well as the nature of the source is: 
T50 — lyr, i.e., the time delay is comparable to the integration time of the experiment, T b s . In 
this case, the limit D6e *C l c holds as the Faraday rotation bound is combined with T50 ~ lyr. 
For a propagated differential energy spectrum dJ oc E^dE, that corresponds to the injected 
spectrum below the GZK cutoff, Eq. (||) tells us that the arrival time distribution is given by 
dJ oc t^-^^dt, for a bursting source with emission timescale T$ <C T50 and a te — E correlation 
with negligible scatter. Note that a continuous source with T$ S> 1 yr would produce a uniform 
distribution of arrival times, notably independent of their energy (see §2.4). 

The observation of the arrival time and energy of two events is sufficient to determine a 
zero-point in time (time of emission), hence the value of t^;, equiv&lently DJ3 rms \fC c . The distance 
can be independently derived by fitting the energy spectrum above the GZK cutoff, provided there 
is enough statistics. Thus, one ends up with an estimate of the distance to the source, the nature 
of the source (e.g., burst vs. continuous emission), and the value of B^ ms l c . Observations of several 
clusters of events would thus "map" the extra-galactic magnetic field. Finally, one can show, from 
T 50 — D61q/2 ~ lyr, that the angular image would not be resolved for a typical detector with 
resolution 50 ^ 0.5°. That typical time delays above tens of EeV may be of the order of a few 
years, is suggested by a detailed likelihood analysis (Sigl, Lemoine, & Olinto 1997) of the three 
pairs of events that were recently reported by the AGASA experiment (Hayashida 1996). 

2.2. Large time delays 

If te S> T b s ~ afewyr, a given source will be seen only on a limited range in energy (Waxman 
& Miralda-Escude 1996); indeed, protons with higher energy have already reached us in the past, 
while protons with lower energy have yet to reach us. These authors have derived the shape of the 
energy spectrum in the limit D6e 3> I a f° r E <J 50EeV where the signal has significant scatter 
AE/E ~ 30%; in the opposite limit DOe *C l c , AE/E ~ 1%. In principle, the measure of the 
width would tell us which limit applies, and would give us some information about £> rms and l c . 
However, this statement depends strongly on the distance, as D9e/I c oc D 3 / 2 . For sufficiently 
large time delays, the angular image of the source could be resolved for E ^ 50 EeV, in which 
case the value of B rms \/Dl c becomes accessible. Moreover, the argument developed by Waxman & 
Miralda-Escude (1996) would allow to place another constraint on these parameters, respectively 
D/l c < 1 or D/l c > 1. 

For sources that are observed above the pion production threshold, the previous statements 
do not apply. However, in this limit, one can use a similar reasoning to obtain an estimate of the 
distance. At a fixed time, only a given range of energies is observed. Intuitively, the larger the 
distance, the more important the pion production, hence the broader the signal. This effect is 
illustrated in Figures [l] and |2[ Figure |l| shows the image of the source in the time-energy plane for 



- 5 - 



a source lying at D = 60Mpc. The correlation te oc E~ 2 is shown as a dotted line. In this plane, 
the effect of pion production is mainly to downscatter from higher energies to lower energies, 
with a trend toward increasing the time delay for a given energy at emission. Figure [2] shows the 
correlation between the width of the signal in energy, as actually seen by the detector, vs. the mean 
energy of the signal. This signal is obtained as a slice in the te — E image, integrated between t 
and t + T b s , where T b s = 5 yr in this case and t is arbitrary; this slice is indicated in Figure [l] by 
the dashed lines. An example of the signal in energy so obtained is shown in Figure || in dashed 
line. The correlation shown in Figure || was obtained by measuring the width of the signal for 
different mean energies, corresponding to different choices of t, and adjusting a straight line fit. 
The shaded areas denote the la range of (numerical) uncertainty. These uncertainties actually 
provide a hint of the actual experimental uncertainties associated with such measurements, as the 
number of particles used in the Monte-Carlo correspond to that expected from a typical UHECR 
source. As Figure [2] reveals, an estimate of the distance 5Mpc <^ D 100 Mpc, could be achieved 
with reasonable accuracy. For D 5 Mpc, the width of the signal becomes dominated by the 
instrumental resolution, whereas for D ^ 100 Mpc, statistics should insufficient. 

2.3. Continously emitting sources 

The case of a continous source, emitting on a timescale Tg, is obtained by folding, on the 
time axis, the angle-time-energy image of a corresponding bursting source, with a top-hat of 
width Tg. In principle, for a given magnetic field configuration, there will be an energy Eg, such 
that Tg = te s - For E 3> Eg, no correlation is expected between arrival times and energies, as 
the arrival time distribution is dominated by the uniform top-hat in this limit. In contrast, for 
E <C Eg the source behaves just like a burst with respect to the observations; since, in general, 
Tg T bs (for radio-galaxy hot spots, for instance), one also has te 3> T b s for E <; Eg. Thus, 
this situation would be analogous to that discussed in §2.2, i.e., the energy spectrum should 
show a cut-off around ~ Eg, as the UHECRs with E <C Eg have not yet reached us, even if 
they were among the first emitted. The transition between these two regimes is only observable 
if lOEeV <J Eg <J a few 100 EeV. Here, the lower bound is given by the requirement that the 
deflection angle in the galactic magnetic field should be less than a few degrees so that observed 
events can be associated with a common source. The upper bound is dictated by the statistics 
required to establish reasonable estimates of the correlations in time and energy. The simulation 
of a possible case is shown in Figure [3|, where the energy spectrum, as seen by the detector, is 
plotted for a total of 50 detected particles. 

The energy spectrum above Eg can be used to estimate the distance D, for a given injection 
spectrum, via a pion production fit (see Figure ||). The angular image should be resolved for 
sufficiently large time delays (see §2.2), and the actual magnitude of te can be derived, as 
te oc D6 e . Therefore, not only can the magnetic field strength B rms y/T c be determined, the 
timescale of emission is also obtained as a by-product. If Tg ~ T t, s , the angular image would not 
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be resolved, but the time delay could be directly measured, as in §2.1, and the above results still 
hold. In the intermediate case, lyr <C T$ <^ 10 3 yr, one could only place an upper limit on the 
magnetic field strength, B xms ^/T c £ 1(T 10 G Mpc 1 / 2 (0 £ ;/O.5 o )(L>/5OMpc)- 1 (^/lOEeV), hence an 
upper limit on te and T$. For reasonable values of the distance, this constraint is already more 
stringent than the Faraday rotation limit. 

If the field strength lies close to the Faraday bound, the range 300(D/30 Mpc) 2 yr <J T$ <^ 
3 x 10 5 (D/30Mpc) 2 yr is within reach. Observation of several sources at various distances would 
enlarge the detectable range of emission time scales and time delays. In case no low-energy cutoff 
is seen down to ~ lOEeV, but the deflection angle is resolved, the above estimates of Ts turn into 
lower limits. The opposite limit is given if Ts <C te for all E up to the value where the experiment 
runs out of statistics. In this case the source behaves like a bursting source and the discussion in 
§2.2 applies. We note that, here, the deflection angle should be measurable, if Ts ^ 10 4 yr, hence 
-BrmsV^c could be measured. 

2.4. Magnetized galactic halo 

In principle, the results of §2.1, §2.2 and §2.3, should be considered as constraints on both the 
EGMF and the halo magnetic field, in the sense that 6 E ~ 9 2 XG + 6 H , where 6xg and Oh denote 
the deflections induced by the EGMF and by the halo field respectively. The halo field has a 
significant effect if its strength is in the range 10~ 8 G <J Bh <S 10~ 6 G and its scale height ^ 1 kpc. 
Present observational constraints are not conclusive (see, e. g., Beck et al. 1996, and Kronberg 
1994); some authors argue for significant halo fields with large scale heights, while others argue 
for smaller scale heights, <; 1 kpc. Therefore, one cannot a priori rule out the possibility that the 
deflection and time delays are dominated by the influence of the halo magnetic field, even if the 
nucleons originate from an extra-galactic source. If the effect of EGMFs are weak, the results of 
§2.1, §2.2, and §2.3 can be applied to the halo magnetic field by substituting the extension of the 
magnetic halo for the distance to the source. The only exception is that, for a bursting source, 
a strong correlation between arrival time and energy is expected even above the GZK cutoff 
irrespective of its distance. This is due to the absence of pion production during propagation 
through the halo when most of the time delay and deflection is accumulated. As well, in case the 
energy spectrum above the GZK cut-off could be observed, for instance for a continuous source 
as in §2.3, or a burst with a small time delay, the absence of pion production would constitute a 
signature of the proximity of this source. 

3. Conclusions 

We have shown that an EGMF of strength 10~ 12 G <J B rms <; 10~ 9 G, and/or a Galactic halo 
magnetic field of strength 10~ 6 G <J Bh <J 10 -6 G, leave distinct signatures in the angle-time- 
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energy images of UHECRs of energy E ^ lOEeV, which could be used, inversely, to probe these 
fields. Obviously, the respective detection or non-detection of a correlation between arrival times 
and energies marks a bursting source with a short time delay te ~ T {, s , or a continuous source 
with emission timescale T$ ^ T b s . A bursting source with a large time delay would be seen with 
no significant correlation between arrival times and energies, however, only in a range of energy of 
width AE/E ^30% (Waxman & Miralda-Escude, 1996). The observation of an energy spectrum 
(integrated over the observational window), up to the highest energies, that shows a cut-off at a 
low energy (still above ~ 10 EeV) constitutes the signature of a continuous source with an activity 
time scale T$ comparable to the typical time delay te at the cut-off energy. Information on the 
actual magnitude of T$ is contained in the high end of the observed spectrum and in the arrival 
directions. The absence of a lower cutoff implies T5 > tiqecV- In the limit of large time delays, 



deflection angles of events around 10 EeV should be measurable, and the value of B rms ^/Dl c could 

be derived therefrom. If both the typical time delay and T$ are smaller than the integration time, 
the whole spectrum above 10 EeV would be "scanned through". In this case, both D and B rms \/J^ 
could be determined. A more quantitative implementation of the effects discussed here should 
involve a likelihood approach. We have performed such an analysis for the three UHECR pairs 
recently suggested by AGASA (Hayashida et al. 1996). Although present data are much to sparse 
to draw any quantitative conclusions, we observed some potentially interesting tendencies (Sigl, 
Lemoine, & Olinto 1997). One of the pairs, for instance, turns out to be inconsistent with a burst 
and comparatively small time delays of the order of a few years may be favored. 

Strictly speaking, the analysis in the present paper applies to models that predict UHECR 
to be nucleons in the relevant energy range. Topological defect models predict a domination of 
7— rays above ~ 50 EeV {e.g., Sigl, Lee, & Coppi 1996). However, due to its electronic content, 
the deflection and delay of an electromagnetic cascade roughly correspond to those of a nucleon 
of the same energy, modified by the relative lifetime fraction of electrons at that energy, typically 
a factor k, 0.5 (Lee 1996). We therefore expect our analysis to reproduce at least the correct 
tendencies also for defect models. 
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Fig. 1. — Contour plot of the UHECR image projected onto the time-energy plane. The distance 
to the source is 60Mpc; B rms = 2 x 1CT 10 G, l c = IMpc, ub = 0. The dotted line indicates the 
energy-time delay correlation te oc E~ 2 as would be obtained in the absence of pion production 
losses. The dashed lines, which are not resolved here, indicate the location (arbitrarily chosen) of 
the observational window, of length T b s = 5 yr. 
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Fig. 2. — Full width at half maximum (FWHM) vs. the average energy, for the observed spectrum, 
above pion production threshold, for large time delays te 3> T ^ s and small deflection D9e *C l c - 
Shown are linear fits to the correlations calculated from time slices through Figure 1 as well as the 
la uncertainty around these fits for different distances D = 5Mpc, D = 30Mpc, and D = 60Mpc, 
as indicated. 



-12- 




E [EeV] 



Fig. 3. — Energy spectra for a continuous source (solid line), and for a burst (dashed line). Both 
spectra are normalized to a total of 50 particles detected. The parameters corresponding to the 
continuous source case are: Tg = 10 4 yr, tioo = 1-3 x 10 3 yr, and the time of observation is 
t = 9 x 10 3 yr, relative to propagation with the speed of light. A low energy cutoff results at the 
energy E$ = 40 EeV where te s = t (see text). The dotted line shows how the spectrum would 
continue if T$ <C 10 4 yr. The case of a bursting source corresponds to a slice of the image in the 
te — E plane, as indicated in Figure 1 by dashed lines. For both spectra, D = 30Mpc, and 7 = 2. 



